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A cooperative network model of sociological interest is examined to determine the sensitivity of 
the global dynamics to having a fraction of the members behaving uncooperatively, that is, being 
in conflict with the majority. We study a condition where in the absence of these uncooperative 
individuals, the contrarians, the control parameter exceeds a critical value and the network is frozen 
in a state of consensus. The network dynamics change with variations in the percentage of contrar¬ 
ians, resulting in a balance between the value of the control parameter and the percentage of those 
in conflict with the majority. We show that the transmission of information from a network B to a 
network A, with a small fraction of lookout members in A who adopt the behavior of B, becomes 
maximal when both networks are assigned the same critical percentage of contrarians. 


The regulatory dynamics of the brain 1], the cardio¬ 
vascular and other physiological systems (2], and indeed 
most biological/sociological networks appear to be poised 
at criticality [3]. The existence of phase transitions is so 
common, in part, because criticality is the most parsi¬ 
monious way for a many-body system, with nonlinear 
interactions to exert self-control. Inhibitory links in neu¬ 
rophysiology and contrarians in sociology are the names 
given to interactions that evoke the disruption of organi¬ 
zation and consensus, thereby suggesting that the well¬ 
being of either the brain or human society requires the 
containment of those negative agents. However, recent 
neurophysiological literature shows that this perspective 
may be overly restrictive, and that a sufficiently large 
concentration of inhibitory links may counter-intuitively 
have the beneficial effect of promoting a ceaseless activity 
g] , a characteristic that can provide evolutionary advan¬ 
tage. 

One of the first explanations of abrupt social tran¬ 
sitions in terms of criticality was made by Callen and 
Shapero in 1974 [5]. They put together the concepts of 
social imitation and critical behavior a generation before 
Gladwell popularized the concept of the tipping point [6]. 
In the sociological phenomenon of interest to us here the 
role of inhibitory links is played by individuals called con¬ 
trarians gj. More recently an interesting connection has 
been established between the action of contrarians and 
the property of frustration found in spin glasses [8J[9]. As 
the term frustration suggests, the action of contrarians 
is found to quench consensus or prevent its occurrence 
in accordance with the sociological conclusions of Croki- 
dakis et al. mu- However, we reach a different conclusion 
and find value in those individuals whose method of de¬ 
cision making are in conflict with the majority. 

Herein the observation made in neurophysiology [J is 
adapted to sociology using the decision making model 
(DMM). The complex network described by the DMM 
implements the echo response hypothesis, which assumes 
that the dynamic properties of a network of identical in¬ 


dividuals is determined by individuals imperfectly copy¬ 
ing the behavior of one another m- The effect of in¬ 
troducing contrarians into a cooperative social network 
is analyzed using a system of coupled two-state master 
equations. Using analytical calculations we show that in 
the presence of contrarians, increased cooperation effort, 
in the form of increased values of the DMM control pa¬ 
rameter, is necessary to achieve consensus. At the same 
time, contrarians may promote a condition of ceaseless 
activity similar to that found in the cognitive context 
g], The maintenance of such a social state of alertness 
supports a kind of flexibility that has survival value. Re¬ 
call the well-known phrase, often attributed to Thomas 
Jefferson, “Eternal vigilance is the price of liberty.” 

With the growing evidence of phase transitions in bi¬ 
ological and sociological systems, recent studies have 
turned to information-theoretic analyses of canonical 
models for more insight. For instance, Refs. EMU 
show that mutual information between elements peaks 
at the phase transition, while Ref. m shows that in the 
Ising model information flow between elements peaks in 
the disordered regime. Following a different but related 
approach, we study the transmission of information from 
one DMM network to another through cross-correlation 
measurements. We demonstrate that information trans¬ 
fer is maximally efficient when both systems are in a crit¬ 
ical condition. This result is in accordance with earlier 
findings mm, although here we point out that a non¬ 
vanishing fraction of contrarians is responsible for bring¬ 
ing the system to the critical point; information transfer 
would be negligible without those individuals in conflict 
with the majority present to disrupt the order in the sys¬ 
tem. 

DMM network dynamics is a member of the Ising uni¬ 
versality class and is found to be useful for describing 
phenomena related to social group behavior (HI 119]. The 
network model is based on the dynamics of single indi¬ 
viduals selecting one of two options. Denoting the two 
options as +1 and —1, the i-th individual generates the 
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stochastic time series = ±1. Using the Gibbs per¬ 

spective this time series is analyzed using the solution to 
the two-state master equation 

+ \g 2 l(t)P 2 \t) 
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where the +1 and —1 options have been labeled 1 and 2 
respectively. The time dependent transition rates g\ 2 (t) 
and (t) determine the production of decision events 
for the z-th individual. At the moment of making a deci¬ 
sion the single individual tosses a coin to decide whether 
to keep the same opinion or to change her mind, hence 
the factors of 1/2 present in the master equation. Al¬ 
though this implies a random decision, the interaction 
with the other individuals may prolong or shorten the 
time necessary to make a decision, thereby generating a 
bias toward one of the two choices. 

For notational simplicity let us describe the behavior 
of the z-th individual omitting the superscript z, while 
keeping in mind for now that this is a single individ¬ 
ual and that there are TV — 1 other individuals in the 
network. The transition rate from state |1) to state |2) 
reads gi 2 = g exp [—K (Mi — M 2 ) / M ], where M is the 
number of individuals linked to the z-th individual, Mi 
is the number of its neighbors in the state |1) and M 2 
the number of its neighbors in the state |2). When the 
interaction coupling parameter K vanishes, the z-th in¬ 
dividual generates a Poisson sequence s(t) with decision 
events generated at the fixed rate g. When K > 0 the 
z-th individual cooperates with her neighbors. That is, 
when the z-th individual is in the state |1) and the ma¬ 
jority of her neighbors share this state, the rate of her 
decision event productions decreases, thereby indicating 
that the z-th individual is likely to remain in the state |1) 
for a more extended time than in the absence of interac¬ 
tion. The same cooperative prescription holds true when 
the z-th individual is in the state |2), leading in this case 
to <721 = 9 exp [-K (. M 2 — Mi) /M] and indicating that 
if the majority of her neighbors are in the state |1) the z- 
th individual makes decisions with a faster rate, thereby 
reducing her sojourn time in the state |2). 

In this article we adopt the all-to-all (ATA) coupling 
condition, which assigns to all the individuals the same 
number of neighbors, M = TV — 1, where TV denotes the 
total number of network members. Since the total num¬ 
ber of members usually satisfies TV ^ 1, we set M = TV. 
Under the ATA coupling condition all the individuals are 
described by only two transition rates. 

The mean field of the network, 
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probability difference as x = pi — p 2 , the master equa¬ 
tion describing the mean field behavior of an ATA DMM, 
consisting of an infinite number of cooperative individu¬ 
als, becomes 

^ = l( e Kx _ e ~K x) _ 9_ {e - Kx + e K x)x (3) 

The equilibrium value of the mean field can be deter¬ 
mined by setting the left-hand side of Eq. © equal to 
zero, which yields the equation for the equilibrium value 
of the mean field: x eq = tanh (Kx eq ). A second-order 
phase transition occurring in the cooperative system at 
K — 1 can be predicted as follows. If we make the as¬ 
sumption that at the phase transition the equilibrium 
value of x is very close to zero, then using the Tay¬ 
lor series expansion of the hyperbolic tangent gives us 
x eq = Kx eq , which is compatible with a small but non¬ 
vanishing solution only for K = 1. 

An individual is a contrarian if she is inclined to 
make a decision that is the opposite of the one made 
by her neighbors [7]. Thus, for instance, the 512 transi¬ 
tion rate would become 512 = g exp [K (Mi — M 2 ) /M], 
and the <721 transition rate would become <721 = 

g exp [-K (Mi — M 2 ) /M\. By following the same line 
of reasoning as that generating Eq. (j3j) we obtain for an 
ATA DMM of contrarians 

f = § (e~ Kv - e Ky ) - §(e*» + e~ Ky )y. (4) 

We use the variable y to denote the mean field of contrari¬ 
ans for the purpose of distinguishing contrarians from co- 
operators. Recall that the variable x is associated with 
individuals that are cooperators. 

When all the people in the network are contrarians, 
the network remains close to the condition of a vanishing 
mean field. In fact, the Taylor series expansion of the hy¬ 
perbolic tangent in the solution now yields y eq = —Ky eq , 
which implies that the network remains fixed at the equi¬ 
librium value y eq = 0, independently of the value K of 
the interaction strength. 

We note however, see Fig. [l] that in a network with 
only a small concentration of contrarians that a phase 
transition occurs with the important symmetry property 
y = —x, indicating that the mean field of contrarians y 
has the same intensity as the mean field of cooperators x, 
but with the opposite sign. We are thus led to examining 
the condition where the global field z, with a fraction 
q of the individuals being contrarians, is expressed by 
z = (1 — q)x + qy , which, using symmetry, becomes z = 
(1 — 2q)x. The master equations for the cooperators can 
then be written 

dx 

— = g sinh [K( 1 — 2q)x\ — gx cosh [K(l — 2q)x\ , (5) 

and the master equation for the contrarians as 


becomes identical to a probability difference in the limit dy 

TV — > 00 where pi = TVi/TV ; z = 1,2. Defining this — = —gsmh[K(l — 2q)x) — gycosh[K(l — 2q)x]. (6) 
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FIG. 1: Mean fields of the ATA DMM networks. The dia¬ 
monds refer to the mean field of cooperators ( x ), the circles 
to the mean field of contrarians ( y ), and the triangles to the 
global mean field (z). N = 10 3 units, g = 0.01, and q = 0.15. 
Note that the mean field of contrarians is opposite to the 
mean field of cooperators. 


Note that the arguments of the exponential functions 
coincide with the global field z, which is perceived by 
both cooperators and contrarians, but contrarians react 
oppositely to that of the majority reaction to the global 
field z. In fact, Eq. §, determining the time evolution 
of y 1 is obtained by replacing K with — K in the expo¬ 
nential functions of Eq. 0 representing the influence of 
the neighbors on the decisions of single individuals. 

We see that, as expected, the equilibrium condition 
generated by Eqs. ([5]) and (pf is x eq = tanh (. K'x eq ) = 
— Ueq , where I\’ = K( 1 — 2q)7 It is evident that a phase 
transition occurs following the same mathematical pre¬ 
scription as in the absence of contrarians with the main 
difference being that the critical value of K is now given 
by 

*•<«> = (7) 

Consequently, the critical control parameter increases in 
value as the fraction of contrarians increases and the in¬ 
teraction effort necessary to make a social decision in 
the presence of 50% contrarians becomes infinitely large. 
Therefore consensus cannot be reached beyond the limit 
of 50% contrarians. 

Fig- i illustrates the formation of a social decision in a 
diverse social network, having a mixture of cooperators 
and contrarians, and the required strength of interaction 
coupling is in agreement with Eq. ([7]). The figure shows 
that when the network is in the supercritical condition 
in the absence of contrarians, the action of an increasing 
number of contrarians has the effect of shifting the net¬ 
work dynamics down towards the critical point. Inverting 



FIG. 2: The mean field of cooperators ( x ) with different frac¬ 
tions of contrarians ( q ) in the ATA DMM network of 10 3 units 
with g = 0.01. Increasing the fraction of contrarians has the 
effect of turning the supercritical into critical condition for a 
convenient fraction of contrarians. 

Eq. 0 to obtain q yields 

=K 1 ~i) • (8) 

indicating that if the interaction strength K is a fixed 
property of the network, there exists a specific fraction 
of contrarians q c that will bring the network to the critical 
point. Thus, the network dynamics can be adjusted to 
operate at the critical point in three distinct ways: 1) 
with no contrarians, a subcritical interaction strength can 
be increased to a critical value; 2) with no contrarians, 
a supercritical interaction strength can be decreased to 
a critical value and 3) for a fixed interaction strength 
above the critical value, the fraction of contrarians can be 
increased to the critical point for the network dynamics. 

West et al. [20] hypothesized that the maximum in¬ 
formation between complex networks occurs when the 
complexity of the two networks, as measured by their re¬ 
spective inverse power law indices, match one another. 
In the present context this hypothesis is tied to the crit¬ 
icality of the network dynamics, since the time intervals 
between decision events display inverse power law behav¬ 
ior at criticality m- Moreover, the interesting discovery 
was made that the transmission of information from a 
driving complex network B to a driven network A be¬ 
comes maximally efficient when the two networks are at 
criticality Ei lig, that is, when they are maximally com¬ 
plex. Here we show that when the percentage of contrar¬ 
ians induces criticality according to the theory presented 
above, the condition of maximal efficiency for the transfer 
of information is realized. 

We consider two identical ATA DMM social networks 
each with N = 10 3 individuals, with their interaction 
strength fixed at K = 1.25. Note, that Eq. ^ predicts 
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FIG. 3: The cross-correlation function ((C))between two pre¬ 
dominantly cooperative networks (the ATA DMM networks 
with g = 0.01 and K = 1.25) versus the fraction of contrari¬ 
ans ( q ). Each network shows a sharp peak in the vicinity of 
the predicted value. 


an interaction strength with 10% contrarians is necessary 
to realize criticality. To connect the driving network B 
to the driven network A 1 we introduce into network A an 
additional 20 “lookout” individuals that track the global 
field of network B m- This choice of coupling between 
the two networks was inspired by a recent experiment 
[22] where the signal from a few electrodes implanted 
in the brain of a rat B is the information transmitted 
directly to the brain of rat A, the whiskers of which, 
quite surprisingly, synchronize with those of rat B. 

Keeping an equal fraction of contrarians q in both net¬ 
works, we evaluate the cross-correlation between network 
A and network B for various q values [21]. Fig. [3] shows 
that the transmission of information from network B to 
network A becomes maximally efficient at q = 0.09. The 
numerical value for the peak of the cross-correlation curve 
is near the theoretical value q c = 0.1, predicted by Eq. Q 
on the basis of the conjecture that criticality maximizes 
the efficiency of the information transport. One possible 
reason for the slight discrepancy between the numerical 
results and the value expected from theory is that Eq. ^ 
does not account for the presence of lookout individuals 
in the DMM network. A fraction of lookout individuals p 
creates the effective interaction strength K. p = K(1 — p) 
[21] . For the numerical simulations we have p = 20/1020 
leading to K p « 1.225. Replacing K by K p in Eq. ^ 
yields q c = 0.092 for the critical fraction of contrarians, 
refining the agreement with the results depicted in Fig. [3] 

We also see from Fig. [3] that when the concentration 
of contrarians tends to vanish, the transmission of infor¬ 
mation between networks becomes very small. This is so 
because the social system falls in the supercritical condi¬ 
tion, which is not resilient and is unsuitable to address 
crucial issues. As a relevant example we have in mind 


the debate on setting a limit to the extraction of energy 
resources so as not to leave the future generations empty 
handed [23]. We believe that criticality corresponds to 
the condition of full democracy that according to the au¬ 
thors of Ref. J23] would be necessary to promote the 
energy sustainability of future generations, and we imag¬ 
ine that this promotion may be realized through a flow 
of information from one social network at criticality to 
another in the same democratic condition. On the other 
hand, a fraction of contrarians larger than the critical 
concentration q c of Eq. ([8]), realizing the subcritical con¬ 
dition, is still compatible with a significant transmission 
of information, due to the distinctly asymmetric shape 
of ( C) as a function of q. Although our theoretical ap¬ 
proach does consider information flow in the same man¬ 
ner as the recent work of Barnett et al. mi, it shows 
that the contrarian-induced disorder favors information 
transmission. 

While the results of the present article lend support 
to the attractive discovery of the neural benefits of in¬ 
hibitory links in cognitive networks [Jj , our findings have 
a distinct sociological significance as well. They suggest 
a kind of equivalence between two apparently quite dif¬ 
ferent forms of complexity, one of neurophysiological in¬ 
terest U and the other of sociological interest, the latter 
pertaining to the ATA DMM used herein. 

The peaking of the cross-correlation function in Fig. [3] 
indicates that the concentration of contrarians within a 
network can be used to establish a form of resonance be¬ 
tween a driven and a driving network, a central result 
of this paper. When the concentration q is assigned the 
critical value forcing the network to transition from a dis¬ 
ordered state to the condition when consensus is possible, 
the two networks establish a kind of synchronization. 

The intensity of the cross-correlation becomes negli¬ 
gible for values of q smaller than the critical value, with 
the network freezing in a rigid consensus with a locked-in 
dependence among individuals, this being a condition of 
flawed democracy [23]. Making q larger than the critical 
value in the limiting case of a very large concentration 
of contrarians has the effect of allowing the individuals 
to recover their Poisson behavior, signaling the statisti¬ 
cal independence between individuals. However, for an 
extended range of concentration of contrarians exceeding 
the critical value, the single units while recovering their 
freedom are still sensitive to criticality and a significant 
amount of information transmission is possible. 

We restricted the discussion of diversity to the mix¬ 
ture of cooperators and contrarians in its simplest form, 
the ATA condition, using a mean field approach. Moving 
away from the ATA condition, to a lattice or generic net¬ 
work, introduces significantly different global behavior. 
Specifically, the role of frustration must be considered, 
as reaching consensus depends also on the topology of 
the distribution of contrarians with some configurations 
favoring and others obstructing criticality. This more 
complicated scenario could be analyzed using a hetero¬ 
geneous mean field theory approach. 
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